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ABSTRACT 

We present accurate new wavenumbers for a set of 13 mid-infrared 
fine-structure lines. The wavenumbers were determined from observations 
of the planetary nebula NGC 7027 and of the red supergiant a Scorpii. 
Most of the new wavenumbers are good to within 0.0025%, or 8 km s _1 . 
We provide details on the measurements and present an analysis of the 
errors. In addition, we present the first observations of hyperfine splitting 
in the [Na IV] 1106 cm" 1 line. 

Subject headings: atomic data — infrared: ISM: lines and bands - 
planetary nebulae: individual (NGC 7027) 



1. Introduction 

The strongest features in the mid-infrared spectra of gaseous nebulae are the 
fine-structure lines of [Ne II], [S III], [S IV], and [Ar III] and the series of unidentified 
PAH or carbon grain features. This set of atomic and solid state features has 
been the center of attention for most mid-infrared galaxy studies to date, and it 
represents the general limit of what can be detected in galaxies from the ground 
with 4-meter class telescopes (e.g. Roche et al. 1991; Kelly et al. 1995). However, 
with the impending launch of the Infrared Space Observatory (ISO) and with the 
development of 8-meter class, infrared-optimized telescopes, new opportunities are 
unfolding in the field of mid-infrared spectroscopy. 

ISO will be capable of observing a large number of lines in bright galaxies, 
including many lines which have never been detected before except in novae. In 
most galaxies, however, the number of detectable lines will be relatively small so 
it will be inefficient to get complete spectral coverage. It is therefore essential to 
have accurate wavenumbers for the lines of interest. In the course of this study, we 
have found that predicted wavenumbers can be in error by as much as several tenths 
of a percent. Thus, in support of the ISO mission and of our own ground-based 
programs, we have undertaken a project to determine accurate new wavenumbers 
for fine-structure lines in the mid-infrared. 
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Our target for most of these observations has been the hot, young planetary 
nebula NGC 7027. This planetary nebula is unusual in that the central star has a 
temperature of roughly 2 x 10 5 K (Kaler & Jacoby 1989), but it is still surrounded 
by a neutral, molecular cloud (see Graham et al. 1993). The kinematics and 
morphologies of the various ions in NGC 7027 could be quite different from each 
other, and the extended, complicated nature of NGC 7027 is the major source of our 
wavenumber uncertainties. 



2. Observations 

We measured a variety of fine structure lines in the mid-infrared spectrum of 
NGC 7027 during observing runs at the NASA Infrared Telescope Facility (IRTF) 
in 1992 August, 1993 June, and 1994 October. The observations were made using 
Irshell, a mid-infrared grating spectrometer (Lacy et al. 1989; Achtermann 1994). 
Irshell has an 11x64 pixel Si: As array, with 11 rows in the spatial direction and 
64 columns in the spectral direction. We used a low dispersion grating in 1992 
August to measure the entire 8-13 /im spectrum of NGC 7027 at a resolution of 
A/AA m 1000. We switched to an echelle grating to measure individual lines during 
the 1993 June and 1994 October runs. The dispersion of the echelle observations 
ranged from 16000-25000 (depending on the grating angle), and with a 2 pixel 
slit width (2" on the 3m IRTF), the resulting resolution ranged from 8000-12500. 
A dome-temperature card, chopped against the sky, was used for flatfielding, 
atmospheric correction, and fluxing (Lacy et al. 1989). For the 1994 October run, we 
installed a test-grade Hughes 20x64 pixel Si: As array. The 20x64 array has slightly 
smaller pixels, and the October 1994 data were measured through a 2 pixel (1.6") 
wide slit at a resolution of 9500-15000. We used the bright star ji Cephei (spectral 
type M2Ia; Hofneit & Jaschek 1982) for flatfielding, atmospheric correction, and 
fluxing of our 1994 October data. 

The ideal approach to this project would be to make spectral maps of each of 
the lines, and we made such a map for the [S III] line in 1993 June (see Figure 1). 
With these data cubes, we could register the spatial and spectral positions of the 
line emission, which would give us the best chance at determining the rest frame 
wavenumber of the emission. As a bonus, we would also get line intensities and 
kinematic information. However, since time was limited, we had to make a number 
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of compromises. First, we observed at only one slit position for each line. Since 
we did not cover the entire source, our information on line strengths was limited 
to one-dimensional line intensity profiles. Second, we did not have the time to be 
precise in our pointing. Our usual practice was to move our N-S slit back and forth 
a few arcseconds E-W across the western emission lobe until we found the maximum 
signal. Our pointing varied because of guiding and tracking errors and because 
the morphologies of the various ions are not identical. Third, we rarely measured 
calibration star spectra. The card calibration is adequate in most cases, but a star 
spectrum is helpful for removing strong telluric features. Fourth, we speeded up 
our observations by opening the slit to two arcseconds even though we could have 
improved our resolution by using a slightly narrower slit. Fifth, we did not make 
long observations. We spent 2-5 minutes on each of the strong lines and no more 
than fifteen minutes on the fainter lines. In Section 3, we discuss the effects of these 
compromises on the accuracy of our wavenumbers. 

The most important result from this paper are the vacuum rest wavenumbers of 
the lines. These values are presented in Table 1, along with estimated la uncertainties 
(see Section 3) and relative intensities. Of particular note are the three hyperfine 
components of [Na IV] at 1106 cm -1 (see Section 6). The reasonable agreement 
between the observational wavenumbers and the wavenumbers predicted by Kaufman 
& Sugar (1986) and Wiese et al. (1966) gives us confidence in our line identifications. 
We present a summary of our wavenumber coverage and non-detections in Table 2. 
Spectra for the lines detected in NGC 7027 are presented in Figure 2. The spectra 
are usually set to zero in regions of very poor atmospheric transmission (<15 %). All 
data were reduced using the Snoopy data reduction program (Achtermann 1992). 



3. Uncertainties in Wavenumber Measurements 

The main sources of error in our wavenumber measurements are summarized 
below. These contributions were summed in quadrature to determine the 
uncertainties listed with the wavenumbers in Table 1. For convenience, we quote 
errors in km s _1 , i.e. as fractions of the wavenumber. 

1. The wavenumber calibration and when possible the dispersion of the mid-infrared 
spectra were determined from the atmospheric lines. We determined the centroids 
of the atmospheric features from our off-position nod frames, and we were able to 
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determine the centroids to within 0.5 - 2 km s 1 . The errors introduced by our 
dispersion solution were always less than 1 km s _1 . 

2. We obtained wavenumbers for the atmospheric lines from the AFGL line list 
(McClatchey et al. 1973), and the wavenumbers are good to well under 1 km s" 1 . 
However, the atmospheric lines are often blended, and their relative strengths depend 
on atmospheric conditions, so there is some uncertainty in the effective wavenumbers 
determined from atmospheric models. When we had to calibrate our wavenumber 
scale from a blend of atmospheric features, our uncertainties ranged up to 4 km s" 1 . 

3. Photon noise led to scatter in the measured wavenumbers of the individual 
emission lines. It is almost impossible to determine the magnitude of this effect 
because of the velocity structure of NGC 7027, but we estimate that the errors are 
negligible for the strong lines and up to 4 km s" 1 for the weaker lines. Since most 
line measurements involve the averaging of detections at several positions along the 
slit, the photon noise is reduced and is expected to contribute no more than 2 km 
s" 1 to the wavenumber uncertainties. 

4. The lines in NGC 7027 do not have simple line profiles. There is often emission at 
more than one velocity from a given position in the nebula. We had limited success 
at deconvolving the line profiles, so we instead determined centroids for the lines. 
We were able to measure line centroids that were accurate to 0.5 - 1 km s" 1 . In some 
cases, bad flatfielding altered the line profiles. In these cases our line centroids could 
be off by as much as 5 km s -1 . The trio of lines at 1106 cm -1 was severely blended, 
leading to 3 - 6 km s" 1 of uncertainty in the line positions. 

5. We hinge our velocity reference frame on the [S III] line. The 3 km s _1 uncertainty 
in the [S III] wavenumber is propagated into the error estimate for all of the other 
lines. The previous best measurements of the [S III] wavenumber were by Baluteau 
et al. (1976), who found a wavenumber of 534.39±0.01 cm -1 based on measurements 
of the Orion Nebula, and by Greenberg, Dyal, & Geballe (1977), who found a 
wavenumber of 534.41 ± 0.03 cm -1 from observations of three planetary nebulae, 
including NGC 7027. We used three methods to estimate the [S III] wavenumber. 
First, by measuring the mean wavenumber over the [S III] map and correcting for 
the heliocentric velocity of 8.8 ± 0.6 km s" 1 given by Schneider et al. (1983), we 
determine an average wavenumber of 534.388 ± 0.004 cm" 1 . Second, by comparing 
the intensity and velocity profiles of the [S III] and [Ne II] emission, we find the best 
match if the [S III] wavenumber is 534.388 ± 0.005 cm" 1 . The rest wavenumber for 
[Ne II] was measured in the laboratory by Yamada, Kanamori, & Hirota (1985), and 
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is 780.424 ± 0.001 cm -1 . Third, by comparing the velocities of the bright ridges 
in our [S III] map to those in the H76o; velocity map by Roelfsema et al. (1991), 
we determine a wavenumber of 534.387 ± 0.005 cm^ 1 . These three values combine 
to form our new [S III] wavenumber estimate of 534.387 ± 0.005 The two [S III] 
slit positions that pass through the center of the nebula have flux-weighted average 
velocity shifts of roughly -1 km s" 1 relative to the nebular mean. 

6. Telescope tracking errors and pointing uncertainties could have a strong influence 
on our line profiles and measured wavenumbers. We see velocity gradients of as 
large as 6 km s _1 per arcsecond along the slit, with a total shift of 15 - 20 km s _1 in 
the line centroids from one end of the slit to the other. In the E-W direction, we 
see a shift of over 10 km s _1 to the red in the flux-weighted mean wavenumber 
along the slit as the slit is moved from west to east across the nebula, with most 
of the shift taking place in the faint outskirts of the nebula. We determined our 
N-S and E-W positions for each line by comparing the intensity and velocity profiles 
along the slit to the [S III] intensity and velocity map presented in Figure 1. These 
comparisons allowed us to determine the slit positions to within 2", which allowed us 
to remove the velocity offsets relative to the [S III] frame with 2 km s _1 uncertainty. 
On the basis of these comparisons, we estimate that the velocity uncertainties due 
to internal motions in NGC 7027 are 2-4 km s _1 . In cases where there was not 
enough flux to determine the profile of the line emission along the slit, we estimate 
our uncertainties to be 4 - 6 km s -1 . 



4. [Fe II] Measurements in Alpha Scorpii 



We observed two [Fe II] lines in the mid-infrared spectrum of the red supergiant 
a Scorpii. These observations will be discussed by Haas et al. (1996). The 
wavenumbers were reduced to the rest frame using the heliocentric velocity of -5.0 
± 0.6 km s _1 determined by Brooke, Lambert, & Barnes (1974) for low-excitation 
metals in the spectrum of a Scorpii. The wavenumber uncertainties are dominated 
by uncertainties in the effective wavenumbers of blended atmospheric lines. The 
[Fe II] wavenumbers are listed in Table 1. 



5. Line Intensities 
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Our slit measurements provide us with intensity measurements for the 
mid-infrared lines in NGC 7027. As can be seen from the [S III] map in Figure 1, the 
brightest emission for the low ionization lines comes from the bright ridges to either 
side of the nebula center (see also Aitken & Roche 1983). Our slit profiles for the 
higher ionization lines are consistent with their having similar spatial distributions. 
Since most of our line measurements were made with our N-S slit crossing one or 
both lobes, we have reasonable estimates for the peak intensities of the lines. These 
peak intensities are listed in Table 1. It should be emphasized that these intensities 
do not all refer to identical spatial positions, and there is a factor of two variation 
in the continuum levels for the various lines. No corrections have been made for 
extinction. These intensities are provided only as a general guideline to the strengths 
of mid-infrared lines in photoionized gas. Line ratios using these numbers could 
easily be off by a factor of two. 



6. Hyperfine Splitting of the [Na IV] Line 



Since the dominant isotope of sodium, Na, has nuclear spin, we suspect 
that the group of lines at 1106 cm -1 are the hyperfine components of the [Na IV] 
3 P2 — 3 P\ transition, which is predicted to lie near this wavenumber. To test 
this hypothesis, we calculated the expected hyperfine splitting and relative line 
intensities. The hyperfine energy splittings are due mainly to the interaction between 
the magnetic dipole moment of the nucleus and the magnetic field at the nucleus due 
to the orbital and spin motion of the electrons. There is a secondary contribution 
due to the interaction between the quadrupole electric field of the nucleus and the 
electrons. Hyperfine splitting has not previously been observed in [ 23 Na IV] but 
was observed by Harvey (1965) in the isoelectronic atom [ 17 I]. To calculate the 
[Na IV] hyperfine splittings, we used the fact that both fine-structure and hyperfine 
structure splittings are proportional to < r~ 3 > to scale from the hyperfine splitting 
measurements of [ 17 I] and the fine-structure splittings of [O I] and [Na IV], 
taking into account the differing nuclear magnetic dipole and electric quadrupole 
moments and the fact that the fine-structure splittings are also proportional to Zj. 
We calculated the expected line intensities using the formulae presented in Townes 
& Schawlow (1955). There are eight hyperfine components in four closely spaced 
groups, from the four hyperfine sublevels of the J=2 fine-structure level. Of these, the 
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two highest wavenumber groups are blended at our resolution, leaving three resolved 
lines. The resulting spectrum (shown as a dotted curve in the [Na IV] spectrum of 
Figure 2) agrees remarkably well with the data, confirming our identification of this 
multiplet of lines. 

We note that 25 Mg, 39 K, and 41 K also have nuclear spin, but the observed 
3 Pq — 3 Pi transitions should have unresolvably small hyperfine splittings. However, 
the hyperfine splitting of the strong 3 P2 — 3 Pi (5.6/im) line of [ 25 Mg V] should be 
resolvable, allowing a determination of the 25 Mg/ 24 Mg abundance ratio. Several 
isotopic K and CI fine-structure lines may also have resolvable hyperfine structure. 

We are very grateful to Matt Richter and Kevin Luhman for their assistance 
with the observations and to the daycrew at the IRTF for their excellent instrument 
support. We thank M. Haas, M. Werner, and E. Becklin for allowing us to present 
the [Fe II] wavenumbers. Special thanks to Jeff Achtermann for developing the 
Snoopy data reduction software and the DSP-based Linus operating system. This 
work has been supported by USAF contract F19628-93-K-0011 and by NSF grant 
AST-9020292. 
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TABLE 1 

Fine-Structure Line Measurements 



Ion 


Transition 


Obs. Date 


Wavenumber 


Peak 


X'ouier 


X,upper 








(cm' 1 ) 


Intensity" 


(eV) 


(eV) 


[Fe II] 


a 4 F 5/2 - a 4 F 7/2 


6-8 Jun 93 


407.84 ± .02 




7.87 


16.18 


[Ne V] 


3 P - 3 Pi 


6 Jun 93 


411.256 ± .006 


1.3(-11) 


97.11 


126.21 


[S III] 


3 Pi- 3 P 2 


4 Jun 93 


534.387 ± .005 


2.6(-12) 


23.33 


34.83 


[Fe II] 


a 4 F 7 /2 - a 4 F 9/2 


6-8 Jun 93 


557.50 ± .02 




7.87 


16.18 


[Mg V] 


3 Pi - 3 Po 


9 Jun 93 


739.581 ± .012 


1.5(-12) 


109.24 


141.27 


[Ar V] 


3 Po - 3 Pi 


20 Aug 92 


763.23 ± .05 


1.3(-12) 


59.81 


75.02 


[Ne II] 


2 P3/2 — 2 Pl/2 


4 Jun 93 


6 780.4238±.001 


2.5(-12) 


21.56 


40.96 


[Na IV] 


3 P2(F=7/2) _ 3 Pi 


10 Jun 93 


1105.88 ± .02 


3 (-13) 


71.64 


98.91 




3 P2(F=5/2) _ 3 Pi 


10 Jun 93 


1106.12 ± .02 


2 (-13) 


71.64 


98.91 




3 p 3 p 
-< 2(F=3/2) — 1 


10 Jun 93 


1106.30 ± .03 


1 (-13) 


71.64 


98.91 




3 P , , , 3 P 
-r 2(F=l/2) _ -fl 


10 Jun 93 






71.64 


98.91 


[Mg VII] 


3 Po - 3 Pi 


10 Jun 93 


1110.00 ± .05 


8 (-14) 


186.51 


224.95 


[Ar III] 


3 P 2 " 3 Pi 


10 Jun 93 


1112.176 ± .015 


5.1(-12) 


27.63 


40.74 


[K VI] 


3 Po - 3 Pi 


10 Jun 93 


1132.52 ± .03 


1.0(-13) 


82.66 


99.89 


[Ar V] 


3 Pl- 3 P 2 


24 Oct 94 


1265.57 ± .04 


1.9(-12) 


59.81 


75.02 


a erg s _1 


cm 2 sr 1 ; not corrected for extinction; uncertainties can be larj 


2;e (see 


Section 5) 



6 wavenumber determined from laboratory measurements by Yamada et al. (1985) 
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TABLE 2 

Non-Detection Summary 



Target Line Wavenumber Coverage Non-Detections 

(cm- 1 ) 



[Fe II] 


406.25 - 


■ 408.60 






[Ne V] 


410.75 - 


• 413.30 


[Ca VI] 


411.5° 


[S III] 


533.8 - 


• 535.1 






[Fe II] 


556.75 - 


• 558.50 






[Mg V] 


735.22 - 


• 741.1 


[Fe III] 


739.1 


[F V] 


742.75 - 


• 746.99 


[FV] 


746.3 


[Ar V] 


759.66 - 


■ 767.12 


[SV] 


767.5" 


[Ne II] 


778.9 - 


• 781.6 






[S III] 


831.70 - 


• 835.18 


[S III] 


832.5 6 


[CI IV] 


848.93 - 


• 853.94 


[CI IV] 


849.6 


[Al VI] 


1095.25 - 


• 1099.78 


[Al VI] 


1097.0 


[Na IV] + 


1104.7 - 


• 1113.0 






[Mg VII] + 










[Ar III] 










[K VI] 


1129.9 - 


• 1133.4 






[Na VI] 


1159.33 - 


■ 1165.20 


[Na VI] 
[Mg IX] 
[CI VI] 


1161.4 

1162.3" 

1165.4 a 


[Ar V] 


1263.4 - 
1268.07 - 


• 1266.0 

• 1270.81 







a line is from an excited term and so is expected to be very weak 
^electric quadrupole transition is expected to be very weak 
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Figure 1: Intensity map for the continuum-subtracted [S III] emission in NGC 
7027. The contour interval is 0.172 erg s -1 cm -2 (cm -1 ) -1 sr -1 . The two contour 
lines in the center of the nebula denote a low flux region. North is roughly 6 degrees 
west of vertical in this plot. Velocities are indicated for each of the points. The 
velocity uncertainties are 4 km s _1 . 

Figure 2: Dereddened spectra for the twelve emission lines detected in NGC 7027. 
Wavenumbers are presented in Table 1 and uncertainties are discussed in Section 3. 
The dotted line in the [Na IV] plot denotes our model fit for the hyperfine splitting. 
The [Na IV] energy splittings were determined by scaling from the hyperfine 
measurements of [ 17 I] by Harvey (1965) and from the fine-structure splittings of 
[O I] and [Na IV]. Relative line intensities were calculated using the formulae of 
Townes & Schawlow (1955). The spectrum has been convolved to have the same 
resolution as the [Na IV] data. 
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